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Abstract 

The anapole form factor of the nucleon is calculated in chiral perturbation theory 
in leading order. To this order, the form factor originates from the pion cloud, 
and is proportional to the non-derivative parity-violating pion-nucleon coupling. 
The momentum dependence of the form factor — and in particular, its radius — is 
completely determined by the pion mass. 
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There has been much recent interest in parity-violating electron scattering on the nu- 
cleon and light nuclei as a source of information about nucleon structure, in particular 
strangeness content. The SAMPLE collaboration has measured the longitudinal electron 
asymmetry for scattering on the proton [|I| and the deuteron at = —q^ = 0.1 GeV^, 
where q is the momentum transferred to the target. Experiments at JLab (HAPPEX 
and GO [§) extend these measurements to higher . 

A class of contributions to the measured electron asymmetry comes from the anapole 
form factor of the nucleon. The anapole is a parity-violating electromagnetic moment 
of a charged particle with spin, related at classical level to the magnetic field induced 
within a torus by a winding wire It vanishes for on-shell photons and thus cannot be 
separated from contact electron-particle operators. Nevertheless, in the case of hadrons at 
low energies, which can be formulated as an effective field theory (EFT), the anapole does 
constitute a class of contributions that is independent of the choice of photon (gauge), 
nucleon and pion fields, and therefore can be examined in itself. The anapole contribution 
to the asymmetry needs to be understood before one can draw definitive conclusions about 
the strangeness current. 

Here we focus on the anapole form factor of the nucleon aX Q <^ Mqcd, where 
Mqcd ~ 1 GeV is the typical mass scale in QCD. In this kinematic regime we can 
perform a systematic expansion of the form factor in powers of Q/Mqcd (times functions 
of Q/rriT^). This expansion receives analytic contributions from short-range physics (rho 
mesons, etc.) and non-analytic contributions from the long-range pion cloud. The latter 
are calculable in a model-independent way in chiral perturbation theory (xPT). We dis- 
cuss the relative order these contributions appear, and show that to lowest order only the 
pion cloud contributes. The anapole form factor at this order is then a calculable func- 
tion of (Q/w^tt)^, completely determined by pion properties and pion-nucleon couplings 
in principle known from other processes. 

The anapole form factor at Q = (the anapole moment) has been calculated before, 
together with subsets of sub-leading contributions 0, ||, and depends on the leading, 
yet poorly-determined, parity-violating pion-nucleon coupling. However, it is not the 
anapole moment per se that is relevant for the above scattering experiments. Because 
it is generated by the pion cloud, the scale that governs the momentum dependence of 
the form factor is set by m^, and not Mqcd- In the chiral limit, the form factor would 
vary wildly between Q = and Q = Mqcd, but even in the real world it could be quite 
different (say, by a factor of a few) at Q ~ 300 MeV (SAMPLE experiment) than at 
Q = 0. For the proper interpretation of the new experimental data one thus needs to 
investigate the momentum dependence of the anapole contribution. 

We present here the leading-order results for the nucleon anapole form factor, which 
turns out to be purely isoscalar. We recover the known result for the anapole moment 
10, P], and further predict its momentum dependence, for which we give an analytic 
expression. This momentum dependence is given by the known pion mass and is an 
example of a low-energy theorem. In particular, the radius of the anapole form factor 
is, as expected, cx and blows up in the chiral limit. These predictions could in the 

future be tested. We will also show that dimensionless factors are such that the scale of 
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the variation is ~ Gm^r, numerically (but not parametrically) close to the rho mass, 
TUp. As a consequence, the form factor at this order does not display much variation with 
momentum. 

Here we denote by iJ^^ parity- violating nucleon current that interacts with the 
electron current —iee-^^e via the photon propagator iDfj,^ = —i{rip_y/q^ + . . .) to produce 
a contribution 

iT = -tee{k')Yeik) D^M)N{p')J:^{q)N{p), (1) 

to the electron-nucleon S matrix. We have = {p — p'Y = < 0- 

We want to evaluate J^^ aX Q <^ Mqcd- Because the nucleon mass is heavy, ttin ^ 
Mqcd, in this regime the nucleon is essentially non-relativistic, and is parametrized by a 
given velocity and spin {S^ = (0, cr/2) in the nucleon rest frame where v'^ = (1, 0)). 
As we are going to see, we can write 

JLiq) = (aoF!^\-Q') + aA'\-q')r3) (S'^q' - S ■ qq^, (2) 

where ^4^(0) = 1 and is the ith Pauli matrix in isospin space, oq (ai) is the isoscalar 
(isovector) anapole moment of the nucleon and aoF^\Q'^) {aiF^\Q'^)) the corresponding 
form factor. 

At Q ~ ttIt,, the resolution of the virtual photon is enough to see pions, and pionic 
contributions have to be taken into account explicitly. Because the delta-nucleon mass 
difference is comparable to the pion mass, the delta isobar is also a relevant low-energy 
degree of freedom. All these contributions can be calculated in a model-independent way 
starting from the most general Lagrangian involving nucleons A^, pions tt and deltas that 
transforms under the symmetries of QCD in the same way as QCD itself. This EFT 
has been described countless times in the literature; here we only highlight the features 
strictly relevant to our calculation, and refer the reader to a good review — such as Ref. 
0— for details. 

Since the symmetries allow an infinite number of interactions, it is imperative to have 
an ordering scheme for the various contributions. Chiral symmetry plays a fundamental 
role here. Such a power counting argument is possible because all strong interactions bring 
in a small scale. Interactions that preserve chiral symmetry involve derivatives of the pion 
field, so they bring to amplitudes powers of the small momentum, or powers of the delta- 
nucleon mass difference; and interactions that break chiral symmetry involve the quark 
masses, so they bring powers of the pion mass. One can then order the strong interactions 
in the Lagrangian, C = J2a according to the chiral index A = d + n/2 — 2, where d is 
the sum of the number of derivatives and powers of the pion mass and of the delta-nucleon 
mass difference, and n is the number of fermion fields ||^. Electromagnetic interactions 
also break chiral symmetry; they do not necessarily contain quark masses, but they are 
proportional to the small charge e. It is convenient to account for factors of e by enlarging 
the definition of d accordingly. It can then be shown that in leading order in the strong 
and electromagnetic interactions, contributions come from 

^%em = liD,^? - + Ntv ■ DN - ^fN{T ■ S ■ D^)N + . . . (3) 
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Here fj, = 93 MeV is the pion decay constant, = {d^ — ieQA^) with Q^^^ = —ie^ab 
and Q*-^-* = (1 + T3)/2, and ". . . " stand for interactions with more pion, nucleon and/or 
delta fields, that are not needed explicitly in the following. Note that at this order the 
nucleon is static and couples only to longitudinal photons. Kinetic corrections and mag- 
netic couplings have relative size 0{Q / Mqcd) and appear in -C^^^g^. The same is true for 
the delta isobar, including the nucleon-delta transition through coupling to a transverse 
photon. The parity-conserving pion-nucleon coupling is given by the axial-vector cou- 
pling of the nucleon; according to naive dimensional analysis Qa = 0(1), which is indeed 

observed, qa = 1-26 (see, e.g., Ref. 0]). A term in >C^2/em provides an 0{{m^/ MqcdY) 
correction that removes the so-called Goldberger-Treiman discrepancy. 

Weak interactions in this EFT have been discussed in Ref. The four-fermion 

interactions between quarks generated by W and Z exchange break chiral symmetry with 
a strength given by the Fermi constant Gp = 1-17 ■ 10~^ GeV~^. The effective operators 
it entails at low energies are proportional to the Fermi constant times the square of a 
mass scale. A natural scale is the pion decay constant, so we assume that these operators 
have coefficients of order of Gpf^ ~ 10~^ times other natural scales 0. There is one 
pion-nucleon interaction with negative index, 

Ci-l = -'^Niry<^hN + ... (4) 

Here h-^^N is the parity- violating pion-nucleon coupling. Again, according to naive di- 
mensional analysis we expect h.„NNfn = GpflMqcD, or h^^j^N = O^GpfirMQCD) ~ 10~^. 
The value of h^^N^ is not well determined (for a review of constraints on parity- violating 
parameters, see, for example, Ref. 0). 

There are various weak interactions with indices 0, 1, and 2. The first pion-nucleon- 
delta coupling, for example, has to contain in the fermion rest frame the 4x2 transition 
matrix Str', in other words, we cannot combine nucleon and delta in a spin zero object. 
Rotational invariance then demands the presence of a gradient to form a scalar, Str ■ V. 
Such interaction can appear only in C^^lak^ and its effects are suppressed by 0{Q / Mqcd) 
relative to those stemming from the Lagrangian (^ . The same reasoning applies to other 
operators. Particularly relevant among higher-order interactions is 

^SL = -\N{a, + h,T,)S,Nd,F^^ + . . . , (5) 
my 

because qq (cii) is a short-range contribution to the isoscalar (isovector) anapole moment. 
From dimensional analysis, di/mj^ = O^eGpf^/m^) = O {eG p / {4:7iy) . Direct short- 
range contributions to the momentum dependence of the anapole form factor first appear 
in d'^laki being further suppressed by 0{{Q/MqcdY)- 

■^There are clearly factors of 47r that might not be accounted for properly this way, but they will not 
affect the relative order of weak interactions. 

'^At lowest order, our parameters relate to those of Ref. [10| (DDH) by hT^NNQA/ f-K — 

rDDH „DDH 
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Figure 1: Diagrams contributing to the nucleon anapole form factor in leading order. 
Solid, dashed and wavy lines represent nucleon, pion, and (virtual) photon, respectively; 
circles and squares stand for interactions from Cf^.^^^ and respectively. For sim- 

plicity only one of two possible orderings are shown here. 

Let us now consider the sizes of specific contributions to J anil)- The first tree level 
contribution is given by the vertex generated by the weak Lagrangian (^. It has a 
size 0(eGi?(5^/(47r)^), and thus contributes 0(e(ji?m^/(47r)^) to the anapole form fac- 
tor (specifically, the anapole moment). The lowest-order one-loop graphs are built out 
of one vertex from the weak Lagrangian (^) and all other vertices from the strong La- 
grangian (^. They are depicted in Fig. [l|. Consider for example the graph in Fig. |l|(a). 
The coupling of the photon to the pion contributes O(eQ), the parity- violating pion- 
nucleon coupling O^GpfT^MqcD), the parity-conserving pion-nucleon coupling 0{Q/ f.,^), 
each pion propagator 0{1/Q^), the nucleon propagator 0{1/Q), and the loop integra- 
tion 0(Q^/(47r)2). We expect the whole graph then to be OieGrMQcoQ/W^), and 
contribute 0{eGpMQCD^'N / {{'^'^)'^Q) to the anapole form factor. One can easily verify 
that Figs. |l|(b),(c) give contributions of the same size. These long-range contributions 
are 0{Mqcd/Q) larger than the most important short-range contribution. The factor 
1/Q represents an infrared enhancement. The anapole moment, which is a constant, is 
O^eGpf^MQCD/^n) ~ 5 X 10"'' e, and diverges in the chiral limit. The anapole form 
factor is OleGpf^MQCD/m^) times a function F{Q'^/ml) = 0(1). 

One can now show that all other contributions (including all delta effects) are of higher 
order. First, other purely short-range contributions start at JO-^lak^ contribute at 
most O^eGpQ'^ /{4:7iMqcd)'^) to the form factor. Second, other long-range contributions 
contribute at most O^eGpQ/ {{4:'7iYMqcd))- Consider other one-loop graphs: they all 
will have at least (i) one vertex from or (ii) one vertex from C^^l^f.- Examples 

are, respectively, graphs where (i) the photon couples to a virtual nucleon or delta mag- 
netically, or causes a nucleon-delta transition; and (ii) the parity- violating pion coupling 
is a nucleon-delta transition. Diagrams containing one insertion of either of these cou- 
plings are formally suppressed by Q/Mqcd compared to the leading-order contributions 
obtained from the charge coupling present in the Lagrangian (^ and the parity- violating 
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pion-nucleon vertex in (^i they first contribute in sub- leading order. A consequence is 
that, although the delta is treated on the same footing as the nucleon, its (virtual) con- 
tributions to the nucleon anapole are all subleading in power counting. Finally, even the 
most important two-loop graphs, made out of CfJ^^^^ and C^^eak^ should be suppressed by 
the usual {Q / in f-!,)"^ associated with loops in xPT. 

The leading contribution to the anapole form factor comes thus from the graphs in 
Fig. |l[ They are evaluated with f ■ g = 0, as a consequence of the static nature of the 
nucleon at this order. In particular. Fig. |l|(c) vanishes. It is straightforward to calculate 
the other diagrams as well. 

The other two sets of graphs combine to give a result in the form of Eq. (Q). We find 
a purely isoscalar result, 

«o = p T- (6) 

ai = 0. (7) 

This result for the anapole moment agrees with Ref. 0. Ref. p contains a partial set 
of sub-leading contributions because it pre-dates the heavy-fermion expansion employed 
here and in Ref. g. Note that Eq. (|) is ~ \/27r/3 larger than our naive estimate. Eqs. 
(H, 0) are predictions of xPT, but unfortunately the value of h.„NN is currently not well 
determined. 

Because the form factor is given by lowest-order loop graphs, it depends on the com- 
bination Q^/m^ only. Because all Q dependence comes from Fig. ^(a), we find that it 
actually is a function of /{2m.,^y: 




This is also a testable prediction of xPT: it is plotted as function of Q in Fig. 

The variation of the form factor with Q is somewhat smaller than expected. This can 
be seen by looking at the anapole square radius (defined in analogy to the charge square 
radius) : 



(r^ ) = -6 

\ an/ ^ 



(0) 



10ml 



(9) 



We see the square radius is smaller than expected by a factor 5. At Q <^ we can write 



^^^(Q^) = 1-^7^ + 0(7^). (10) 

This radius approximation to the form factor is also displayed in Fig. 

The accidentally small square radius makes the mass scale that governs the Q depen- 
dence near (5 = larger than 2771-,^. Data on the proton Sachs form factors are well fitted 



'The Lagrangian of Ref. [pi can be obtained from ours by a redefinition of the sign of the pion field. 
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Anapole Form Factor 




Figure 2: The isoscalar anapole form factor F\ as function of Q: xPT in leading or- 
der, Eq. (H) (solid line); quadratic approximation, Eq. ( p!Ol) (dotted line); and dipole 
approximation, Eq. (|11|) (dashed line). 

by a dipole profile with a mass scale close to m^. A dipole form factor with the correct 
anapole square radius is 

with the mass scale M = 2\/T0mT^ = 880 MeV, instead of 2711^^. M is numerically close to 
rrip, but this is clearly a coincidence that would be destroyed if the quark masses where 
much smaller than what they are. In fact, short-range contributions such as that from 
the rho appear only in the counterterms at higher orders, so M has no obvious connection 
to rUp. The dipole approximation is also shown in Fig. 0. It improves over the radius 
approximation but clearly Eq. (H) is softer than a dipole. 

In conclusion, we have shown results for the anapole form factor of the nucleon in 
leading order in xPT. We have not yet examined extensively the form factor in next 
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order. There undetermined short-range isoscalar and isovector parameters appear, and 
the anapole moment can no longer be predicted in a model-independent way. However, 
the momentum dependence continues to be determined by the pion cloud — with nucleons 
and deltas in internal lines — and an improved low-energy theorem can be derived. Under 
the assumption that higher-order results are afflicted by the same dimensionless factors 
seen above, the error in Fig. |^ at momentum Q would be ~ Q/nip. Within such an 
error, the approximation of using the anapole moment instead of the full form factor 
would be good to about 15% in the SAMPLE experiment, but larger changes appear 
at momenta relevant to the JLab experiments. On the other hand, if in next order the 
dimensionless factors that affect the variation turn out closer to 1 than to 1/5, the 
momentum dependence from sub-leading order could be (accidentally) comparable to the 
leading order considered above. Only an explicit calculation can verify if the contributions 
that are formally subleading are indeed numerically small. We are currently investigating 
these issues [O . 
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Note added 

After this work was completed, it was pointed out to us that the nucleon form factor in 
leading order had previously been calculated in |T2|. 
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